L Lyapunov Exponent of Quadratic Function Source Code 



> restart; 

Calculate the Lyapunov exponent for f(x) = r*x*(l-x) 

Use r = 3.80, 3.81, 3.82, 3.83, 3.84, 3.85, 3.86, 3.87, 3.88, 3.89. 

Use the initial condition for intial APACHE III probability of mortality (example 50%) xO 
:=0.5. 

The variable 

s = s(x0,i) := ln(|df(xi)|) + s(i-l) as the sum of the logarithms of the first i iterates. 
We take s(x0,0) = 0, so when we take xO = 0.5 we do not calculate the derivative at 0.5. 
The variable 

h = h(x0,i) := s(x0,i)/(i+l) 
is the estimate of the Lyapunov exponent starting with the initial condition xO after i 
iterates. 

(In the text, we write the dependence of s and h on the initial condition xO and the iterate 
i> 

but in the calculation we only use the variables s and h.) 

For the first n 1=20000 iterates the calculation for h(x0,i) is not printed out. 
For the next n2 = 20 iterates the calculation for h(x0,i) is printed out, which is the 
estimate for the 
Lyapunov exponent. 

f The function f(x) = r*x*(l-x) 

df The derivative df(x) = f (x) 

r: Parameter value. 

nl The number of iterates for which the estimate of h(x0) is not displayed 

n2 The number of iterates for which the estimate of h(x0) is displayed 

xO: The initial condition 

> r : = 3 . 81; 

f : = x -> r*x* (1-x) ; 

df := x -> r - 2*r*x; 

nl := 2 0000; 

n2 := 100; 

>x0 := 0.5: 

x : = xO: 

s := 0: 

for i from 1 to nl do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 

od: 

printfC\n\t r = %1.4f \n\n**, r) ; 
printfC\t xO = %1.4f \n\n"\ xO) ; 
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printfTXt i \t\t h(xO,i) \n^) ; 
for j from 1 to n2 do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 
h : = s/ ( j +nl) : 

printfC\t %d \t %1.4f \rT , j+nl, h ) 

od: 

>r := 3.81; 

f : = x -> r*x* (1-x) ; 

df := x -> r - 2*r*x; 

nl := 20000; 

n2 := 100; 

>x0 : = 0.2: 

x : = xO: 

s : = 0 : 

for i from 1 to nl do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 

od: 

printfr\n\t x0 = %1.4f \n\n^, xO) ; 
printfT\t i \t\t h(x0,i) \n^); 
for j from 1 to n2 do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 
h := s/ ( j+nl) : 

printfT\t %d \t %1.4f , j+nl # h ) 

od: 

>r := 3.81; 

f : = x -> r*x* (1-x) ; 

df := x -> r - 2*r*x; 

nl := 20000; 

n2 := 100; 

>x0 := 0.3: 

x := xO: 

s := 0: 

for i from 1 to nl do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 

od: 

printf(^\n\t xO = %1.4f \n\n\ xO) ; 
printfC\t i \t\t h(x0 / i) \n^); 
for j from 1 to n2 do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 
h := s/ (j+nl) : 

printfT\t %d \t %1.4f \n\ j+nl # h ) 

od: 
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V. * 



>r := 3.81; 

f := x - > r*x* (1-x) ; 

df := x -> r - 2*r*x; 

nl := 20000; 

n2 := 100; 

>x0 := 0.49: 

x := xO: 

s := 0: 

for i from 1 to nl do 
x := f (x) ; 

s := s + In (abs (df (x) ) ) : 

od: 

printfr\n\t xO = %1.4f \n\n^, xO) ; 
printfp\t i \t\t h(x0,i) \n"); 
for j from 1 to n2 do 
x := f(x); 

s := s + In (abs (df (x) ) ) : 
h : = s/ ( j+nl) : 

printfTXt %d \t %1.4f \n\ j+nl, h ); 

od: 

> 
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